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ABSTRACT 

Extrapolation  techniques  including  the  geos trophic  and  gradient 
thermal  wind,  average  wind  shear  from  United  States  Weather  Bureau 
msxi.mum-wind  analyses,  and  persistence  were  statistically  compared 
for  January  1°60  over  the  united  States.      Calculations  were  made  for 
layers  contained  between  200  mb  and  1+0,000  or  15,000  feet,   the  thick- 
nesses averaging  2100  or  6^00  feet,   respectively.     Graphic  results  are 
presented  for  both  speed  and  direction  errors  for  each  type  of  extrapol- 
ation. 

It  is  shown  that  low  wind  speed,  cross-isobar  flow,  and  large 
curvature  are  associated  with  significant  extrapolation  errors. 

The  author  is  grateful  for  the  patience  and  expert  counsel  of 
Professor  it.  J.   Renard  of  the  United  States  Waval  irostgraduate  School, 
Monterey,    California. 
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1.   Introduction 

The  progress  of  jet  aircraft  to  higher  and  higher  operat- 
ing altitudes  has  created  a  demand  for  more  precise  upper- 
atmosphere  information.   At  present,  there  is  a  particular 
need  for  an  accurate  description  of  the  wind  field  in  and 
about  the  tropopause.   The  meteorological  parameters  of  most 
concern  to  aircraft  operating  at  these  high  altitudes  are  the 
enroute  winds  and  clear  air  turbulence.  The  former  is  normally 
obtained  from  soundings  or  other  types  of  direct  wind  (or 
pressure)  observations. 

Actual  upper  level  wind  reports  are  frequently  missing 
for  such  areas  as  oceans,  unpopulated  continental  sectors,  or 
enemy  held  territory  in  wartime.   This  deficiency  may  be 
corrected  by  vertical  extrapolation  of  wind  from  a  level  of 
known  data  by  employing  various  aids  such  as  the  horizontal 
temperature  field,  empirical  values  of  the  wind  shear,  or 
layer  of  maximum  wind  analyses,  to  i  atne  a  Eew   Although  each 
of  these  parameters  have  been  studied  individually  by  others, 
no  comprehensive  comparison  of  vertical  extrapolations  of  the 
horizontal  wind  at  high  altitudes  has  been  published. 

This  shortcoming  prompted  a  comparative  analysis  of  the 
following  vertical  extrapolation  techniques  in  the  vicinity 
of  the  tropopause:   geostrophic  thermal  wind,  gradient  thermal 
wind,  average  vertical  shear  of  the  actual  wind,  and  persistence 


2.   I'revious  extrapolation  studies 

Various  methods  are  presently  employed  to  determine  the 
horizontal  wind  field  above  the  highest  reported  wind.   The 
simplest  extrapolation  procedure  considers  the  wind  vector  to 
remain  constant  with  height.   This  is  the  best  technique  to 
use  for  short  vertical  distances  (500  feet)  or  in  areas  where 
little  change  in  the  wind  vector  with  height  is  indicated. 
Determination  of  the  latter  areas  is  dependent  upon  an 
accurate  subjective  analysis  of  the  wind  field  below  the 
desired  level,  on  the  dynamics  of  the  wind  field  as  associated 
with  the  temperature  field  and  the  location  of  the  jet  axis. 
The  condition  of  Z\ V/Zs&  =   O        is  least  likely  to  be  met 
a  few  thousand  feet  above  and  below  the  tropopause  and  level 
of  maximum  wind. 

The  geostrophic  thermal  wind  is  often  used  for  extra- 
polation thrpugh  layers  in  which  the  wind  field  is  changing 
rapidly.   This  method  requires  the  determination  of  the  mean 
horizontal  temperature  field  for  the  layer.  Lacking  a  mean 
temperature  field,  the  isotherms  at  some  level  within  the  layer 
are  often  substituted  for  that  of  the  former. 

A  few  synoptic  studies  of  the  geostrophic  thermal  wind 
have  been  made.   In  a  recent  survey  by  Cunningham  (_ 2 J  accurate 
wind  and  temperature  data  were  obtained  from  a  specially 
instrumented  B-29  flying  at  280  and  383  mb  over  the  United 
States  on  tracks  perpendicular  to  the  jet  stream.   Percent 
departures  of  the  geostrophic  and  gradient  thermal  wind 


extrapolations  from  the  observed  shear  decrease  exponen- 
tially     increasing  distance  (d)  over  which  the  horizontal 
temperature  gradient  was  measured.  Results  were  only  slightly 
better  for  the  gradient  computations.  Cunningham  considered  the 
lack  of  a  large  difference  between  the  two  methods  due  to  the 
small  isobar  curvatures  in  the  sample  (maximum  curvature 
approximately  1/1000  nautical  miles). 

To  illustrate  the  advantage  of  using  a  mean  wind, 
Cunningham  averaged  the  wind  over  d  before  vectorial ly  adding 
the  geostrophic  thermal  wind  to  obtain  the  upper  level  wind. 
The  departure  of  this  wind  from  the  actual  wind  varied  inversely 
with  d,  and  decreased  from  11  percent  at  33  miles  to  2.7 
percent  at  181  miles.  Even  better  results  (two  percent  improve- 
ment) were  achieved  when  using  a  mean  verification  wind  as 
averaged  over  d,  rather  than  a  spot  wind. 

In  order  to  understand  the  wind  distribution  about  the 
jet  core  and  perhaps  to  subjectively  adjust  vertical  wind  shear 
extrapolations,  a  knowledge  of  the  jet  stream  is  essential. 
Several  empirical  investigations  of  vertical  shear  in  the 
vicinity  of  the  jet  stream  have  been  recently  published. 
Endlich  (_3J  has  developed  a  model  cross  section  through  the 
jet  stream  for  use  as  an  aid  to  forecasters  in  filling  in  areas 
of  missing  data.  It  provides  an  average  wind  distribution  in 
a  plane  perpendicular  to  the  jet  in  terms  of  percent  of  the 
speed  at  the  jet  core. 


In  a  report  on  the  layer  of  maximum  wind,  Reiter  j_6j 
indicates  that,  among  other  things,  there  is  a  lack  of  systema- 
tic  dependence  of  the  vertical  decay  of  the  wind  on  the 
magnitude  of  the  maximum  wind.   However,  he  does  show  percent 
vertical  decay  as  a  function  of  the  horizontal  distance  from 
the  jet  core.   A  later  United  States  Navy  AROWA  study  J  10; 
obtained  jet  stream  vertical  shears  that  showed  an  increase 
with  jet  core  speed.  The  magnitude  of  vertical  shears  above 
the  core  were  indicated  to  be  as  much  as  30  percent  greater 
than  those  below.   However,  a  United  States  Air  Force  study 
[  9j  of  winds  over  100  knots  indicates  the  maximum  vertical 
shears  to  lie  below  the  maximum  wind  level. 

A  particularly  useful  description  of  the  jet  streams  over 
the  Pacific  Ocean  is  provided  by  Serebreny  [8j  who  presents 
numerous  rules  of  thumb  from  statistical  studies  of  these 
jets.  A  more  comprehensive  theoretical  presentation  of  jet 
stream  information  is  provided  by  Riehl  et  al  [7j  . 


3 .   Theoret i 

The  equation  of  horizontal  motion  in  natural  coordinates  can 
be  written  as  |_5J 

-V«2  i  (vi ,      )y  *       %■■  *L  (') 

This  becomes  the  gradient  wind  equation  for  V  =  0,  and  as  a 
special  case,  the  geostrophic  wind  equation  when,  additionally, 
Kh  =  0. 

In  a  manner  similar  to  the  development  by  Forsythe  '4'  , 
equation  (1)  is  differentiated  with  respect  to  p  to  obtain 


Interchanging  the  Vp   and  ?L_   operators,  equation  (2)  becomes 

Substituting  the  hydrostatic  equation  (   - ■  ■ >  w     '  —    j        in 
equation  (3)  and  dividing  by  < 

£    ^-VVIViH-lyiu^iv-xt      (4) 


-r  \ 


Taking  the  logarithm  of  the  equation  of  state  (p<*  ~ 

In  p  -r  |k  <*    In  'v  •+"  I"  ;'  i  (5") 

Applying  the  operator  "V^   to  equation  (5)  yields 

;    5    -     =  -~    V.   i 


tituti  .  luation   (4)  a 

"^+(WM  .-*(*&)=  -3       K     (7) 

utial  acceleration,  V,   in  the  first      ,,f  (7)  is 
;i  to  the  component  of  the  pressure  force  in  the  eirection 
of  the  wind  and  is  the  most  difficult  to  evaluate.   The  first 
ay  be  expanded  by  sett;  . 

w      .'    is  the  angle  between  the  geostrophic  ana  actual 
Lions.    Thus  equation  (7)  becomes 

?  #      2  2  <?  2    * x    .-■'  *  - 

accelerations  exist,  A     equals  20  to  30 

3  i 

-—    it    frequently  three      es  per  the 

occ.  as  nine  degrees  per  thous,        |.,J  . 

Lehl  (_  7J  have  found  A      to 

cent  of  the  time  and  greater  than  20 

>nlj  .   percent  of  the  time  in  the  vicinity  of  the         am. 

so  rmination  of  ,3    (much  less  ?!/2      )  from  routine 

/  ol 

uour  ana    is  is  questionable.   However,  when  values  of^' 
definitely  greater  than  zero  appeared,  they  were  tabulated  as 
aii  indication  of  large  V. 


this  study,  201        been  selected  as  the  known  level. 
Little  moisture  exists  at  this  height  and  the  temperature  averages 
approximately  -55°C  (218°K)  over  the  United  States  for  the  period 
selecteu  for  study.   The  temperature  field  is  assumed  invariant 
from  200  mb  to  the  extrapolated  level  with   V„  lv  -  Vp  T 
(temperature  measured  at  200  mb ) .   With  the  further  assumption 
that  Kjj  and  V  are  identically  zero,  equation  (7)  becomes  the 
geostrophic  thermal  wind  equation 

expressed  in  its  scalar  working  form  as 

1     <j      ^°    T       /ooo  i  \       ' 

where  V^  is  obtained^  by  applying  a  geostrophic  wind  scale  for  a 
contour  interval  of  200  feet,  to  isotherms  at  a  5°C  interval. 

The  gradient  thermal  wind  equation  is  determined  from 
equation  (7)  by  setting  V  equal  to  zero  and  using  the  above 
definition  of  the  vector      ■■■  i*    4     The  equation  is  shown 
here  in  a  finite  difference  form  as  a  change  through  a  layer 
/Cs    z,  where  K,  is  assumed  constant. 

(a)  (eO  .  ,   (c) 

4-  f 


Equation  (12)  is  expressed  in  this  form  for  ease  in  handling 

the  term  containing  the  geostrophic  thermal  wind  vector.   Figures 

1  and  2  illustrate  schematic  wind  extrapolations  using  equations 
(10)  and  (12),  respectively. 

The  horizontal  trajectory  curvature,  Kh,  is  obtained  from 
Blaton's  equation 

Kh-  K((i-^-»)  M 


^  fA) 


Vj  =  ^i  *  ^v^ 


Figure  1.   Schematic  Geostrophic  Thermal  Wind  Vector  Diagram 


A\/^  (A) 


Vj.  =  V^  +  AV^r  =  ZWgi  irW  +  iti 


Figure  2.   Schematic  Gradient  Thermal  Wind  Vector  Diagram 


4.   Computation  procedures 

All  data  wei         .  .  from  1200Z  January  I960  United  States 
Weather  Bureau  facsimile  charts  and  concurrent  upper  wind  reports 
200f nb  was  selected  as  the  level  of  known  data  because,  averaging 
39,000  feet  for  this  data  sample,  it  was  considered  to  be  suf- 
ficiently close  to  both  the  level  of  maximum  wind  and  the  present 
jet  aircraft  maximum  operating  altitude.   At  each  station  extra- 
polations of  winds  to  40,000  and  45,000  feet  were  compared  to 
reported  winds  at  those  levels.   United  States  Weather  Bureau 
tropopause-height  and  maximum-wind  analyses  for  this  period  were 
also  utilized. 

a.   Thermal  wuid  extrapolations 

Stations  were  chosen  where  a  well  defined  thermal  gradient 

existed  from  which  the  wind  factor,  V*t>  could  be  measured  with  a 

geostrophic  wind  scale  as  described  i;.  section  3.   V't  is  inulti- 

i 

plied  by  .llZ\z  per  thousand  feet  to  obtain  (A)  in  equation  (12) 
and  Figure  1.   V^  is  the  vector  wind  at  200  mb .   The  wind  at 
the  upper  level  (V^_)  is  obtained  by  the  vector  (graphical) 
addition  of  Z^  Vgs  to  V]_ . 

For  term  (13),  isobar  curvature  (Ks)  is  measured  by  the  use 
of  overlay  arcs  corrected  for  map  scale  variation.   The  200-mb 
wind  (V;l)  and  associated  pressure-system  movement  (c)  are 
determined  and  entered  in  equation  (13)  to  yield  trajectory 
curvature,  K^>  considered  constant  through  the  layer. 


The  scalar  Vgr   is  approximately  equal  tc  V,  the  actual 
mean  wind  of  the  layer.   As  a  first  estimate  of  V  ,  an 

10 


approximate  vectoi  mean  of  V&  and  V^  may  be  obtained.   An  improve- 
ment is  provided  by  the  vector  mean  of  Vr  and  Vi  . 

Term  (A),  /\Vgs,  is  usee  for  the  first  estimate  of   AVgr. 
The  coriolis  parameter  (f)  is  taken  from  published  tables  ["II  . 
Although  (B)  is  a  vector  parallel  to  /\Vgr,  it  may  be  drawn 
parallel  to   AVgs  as  a  first  approximation. 

In  term  (C)    /\Vgr  is  initially  considered  equal  to  the 
scalar  magnitude  of   /\Vgs.   Vector  (C)  parallels  Vgr ,  approxi- 
mately by  V,  which  may  be  obtained  as  previously  described  or  by 
trial  and  error. 

The  resultant   AVt-r  1S  reentered  in  (B)  and  (C),  Vgr  is 
recomputed  and  a  second  approximation  to   /^Vgr  calculated.   This 
is  generally  sufficient  to  produce  the  desired  accuracy  in 
A  Vgr. 

b.  s.      Lnd  extrapolations 

Over  the  United  States,  the  facsimile  maximum-wind  anal- 
yses provide  the  height  and  velocity  of  the  maximum  wind  and 
the  vertical  wind  shear  in  the  layer  containing  the  level  of 
maximum  wind.   (This  level  is  hereafter  referred  to  as  the 
LMW).   These  shears  are  a  subjectively  determined  average  of 
the  plotted  shears  10,000  feet  above  and  below  the  LMW. 

The  shear  obtained  from  the  chart  is  multiplied  by  the 
thickness  of  the  extrapolation  layer  (in  thousands  of  feet).   This 
gives  /\Vs  which  is  added  to  or  subtracted  from  the  wind  at  the 
known  level  (200  nib),  depending  upon  whether  it  is  below  or 
above  the  LMW.   For  extrapolations  crossing  the  LMW,  a  net   /\  Vs 


11 


is  obtained  from  the  difference  between     Vs  from  200  mb  to 
the  LMW  and  t    .  Vs  from  the  LMW  to  the  desired  upper  level.   The 
wind  direction  is  considered  to  remain  constant  through  the 
layer. 

A  shear— wind  computation  was  made  for  all  stations  for 
which  a  thermal  wind  computation  was  made. 
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5.  Graphic  analysis       putations 

The  objective  of  this  study  was  to  statistically  consider 
the  error  distribution  of  various  methods  of  vertical  extra- 
polation of  the  horizontal  wind  field  and  its  relation  to 
other  parameters,  and  thus  determine  the  most  accurate  extra- 
polation procedure.  A  graphic  representation  was  selected  as 
best  meeting  these  requirements. 

a,   Error  discussion 
Error  as  used  herein  is  defined  as: 

(1)  tiie  angle  between  the  extrapolated  wind  direction  and 
the  actual  wind  direction  (from  teletype  reports)  at  a  given 
level.  A  plus  error  means  that  the  extrapolated  wind  had 
turned  through  a  greater  angle  (and  in  the  same  direction)  than 
the  actual  wind  for  the  given  layer,  the  bottom  of  which  is  at 
200  mb .  A  minus  error  represents  too  little  turning  (or  turning 
in  the  wrong  direction)  of  the  extrapolated  wind.   Therefore, 
for  extrapolations  using  persistence,  all  errors  are  minus. 
However,  since  the  distributions  of  the  thermal^ wind-^direction 
extrapolations  are  nearly  normal,  a  clearer  comparison  is 
obtained  by  distributing  the  persistence  extrapolation  errors 
symmetrically  about  an  error  of  zero,  i.  e.  for  any  given  error, 
half  the  cases  are  considered  plus,  half  are  taken  as  minus. 

(2)  the  difference  between  the  extrapolated  wind  speed  and  the 
actual  wind  speed  (from  teletype  reports)  divided  by  the  actual 
wind  speed  at  the  level.   Thus,  the  result  is  expressed  in 
percent  of  the  actual  wind  speed. 
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It  must  be  realized  that,  at  times,  extreme  errors  may 
arise  solely  due  to  errors  in  observing  and  reporting  the  winds. 
Reite.      ;as  shown  how  the  method  of  determining  upper  winds 
can  produce  the  unrealistic  oscillations  sometimes  found  in 
vertical  wind  profiles.   He  also  described  how  the  abbreviated 
teletype  codes  can  remarkably  alter  the  actual  upper-level 
sounding  curve.   His  illustrations  particularly  emphasize  the 
desirability  of  averaging  and  smoothing  to  obtain  representative 
wind  profiles. 

Another  error  arises  when  the  neglected  term,  the 
tangential  acceleration  becomes  large.  Subject  to  errors  in 
analysis,  cross~isobar  flow  is  an  indication  that  this  accelera- 
tion is  significant.  To  determine  the  effect  of  this  term, 
the  number  of  clearly  apparent  cross-isobar  angles  greater  than 
15  degrees  have  been  recorded  in  boxes  (  LJ  )  on  the  graphs  which 
follow. 

The  number  of  winds  under  20  knots  have  also  been  recorded 
in  circles  (  £)  )  on  the  graphs.  At  this  low  wind  speed, 
large  errors  in  direction  and  percent  speed  result  from  small 
vector  errors. 

b.  Data 

Information  from  150  soundings  were  utilized  to  obtain 
107  layers  whose  thickness  averaged  2100  feet  (range  from 
900  to  3300  feet)  and  149  layers  whose  thickness  averaged  6500 
feet  (range  from  3300  to  8400  feet).  Two  thermal  wind  extra- 
polations and  a  shear  extrapolation  were  compared  with  persis- 
tence for  each  layer,  considering  direction  and  speed  separately. 
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Of  the  150  reports,  21  showed  a  near-zero  the;      adient 
at  the  200-mb  level,  45  had  trajectories  which  were  nearly 
straight  (curvature  less  than  1/2500  nautical  miles),  22  cases 
were  associated  with  anticyclonic  curvature,  and  62  with 
cyclonic  curvature. 

The  following  factors  were  considered  for  their  effect  on 
the  extrapolation  error: 

1.  Cross -isobar  angle 

2.  Location  of  the  jet  maximum  in  relation  to  the  station 
(horizontal  distance  and  direction) 

3.  Tropopause  height 

4.  Level  of  maximum  wind 

5.  Whether  wind  speed  was  increasing  or  decreasing  with 
height 

6.  Angle  between  the  thermal  wind  and  actual  wind 

7.  Magnitude  of  wind  speed 

8.  Magnitude  and  algebraic  sign  of  trajectory  curvature 

9.  Magnitude  of  the  thermal  gradient 

The  graphs  that  follow  represent  the  most  significant  results 
of  the  above  considerations. 
Figure  3. 

This  graph  presents  frequency  of  errors  in  the  extra~ 
polated  angular  turning  of  the  wind  for  layers  averaging  2100 
feet.  Little  deviation  from  persistence  is  shown  by  the  geo- 
s trophic  or  gradient  thermal  wind  errors,  since  87  percent  of 
the  geostrophic  or  gradient  thermal  wind  computations  lie  within 
-  15  degrees  of  the  actual  wind.   Of  the  persistence  cases, 
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84  percent  fall  within  this  range. 

The  slight  skew  toward  negative  angle  error  indicates  that 
the  computations,  on  the  average,  underestimate  the  actual 
angular  turning.  A  further  analysis  of  these  winds  indicates 
that  this  slight  skewness  is  due  to  cases  in  which  the  wind 
speed  decreases  with  increasing  height. 

Considering  only  the  cases  of  K^ > 1/2000  nautical  miles, 
where  C±   Vgs  exceeds  two  knots  per  thousand  feet,  does  not 
change  the  substance  of  the  foregoing  conclusions. 

As  would  be  expected  for  extreme  angular  errors 
(>  I  40  degreesj ),  there  are  relatively  large  percentages  of 
low-wind-speed,  cross-isobar  and  large-curvature  cases. 

Figure  4 

This  figure  presents  the  same  type  of  information  displayed 

in  Figure  3,  but  only  for  extrapolations  through  layers  aver*g- 

+  + 

ing  6500  feet.   For  -  15  degrees  and  -  25  degrees  errors, 

respectively,  the  included  percentage  frequencies  are:   63  and 
76  percent  for  persistence,  70  and  84  percent  for  geostrophic 
and  73  and  88  percent  for  gradient  thermal  wind  extrapola- 
tions.  The  latter  represents  a  significant  improvement  over 
persistence,  although  slighly  skewed. 
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The  severest  skewing  is  observed  for  the  geostrophic  cases, 
in  contrast  to  figure  3„   However,  again,  this  effect  i 
the  cases  of  winds  which  show  decrease  with  increasing  height. 
At  the  extremes  of  the  distribution,  low  wind  speeds  account  for 
a  large  portion  of  the  error. 

Large  curvature  is  present  in  most  of  the  cases  associated 
with  errors  greater  than  20  degrees.  Note  that  the  mode  for 
cases  of  large  curvature  lies  at  -10  degrees.   Even  so,  74 
percent  of  the  geostrophic- and  82  percent  of  the  gradient* thermal 
wind  cases  in  this  category  lie  within  -  25  degrees  error. 
Figure  5 

The  distribution  of  the  percent  error  in  the  extrapolated 
wind  speed  from  the  actual  (reported)  wind  speed  is  shown  here, 
for  layers  averaging  2100  feet.  All  four  types  of  computations 

peak  together  at  zero  with  little  skew  except  possibly  the  A  V$T 

+  f- 

wind  cases.  Within  -  15  percent  and  -  25  percent,  respectively 

lies  66  and  80  percent  of  the  persistence  cases,  68  and  85  pe 

cent  of  the  shear  cases,  71  and  89  percent  of  the  geostrophi* 

and  75  and  90  percent  of  the  gradient-thermal  wind  cases. 

Persistence  and  shear  show  the  maximum  large  error  samples. 
As  expected,  they  are  associated  with  a  correspondingly  large 
number  of  cross -isobar  and  low-wind-speed  cases.   Cases  with 
large  curvature  are  considered  on  the  next  two  graphs. 
Figure  6 

This  figure  depicts  the  error  distributions  for  those 
extrapolations  in  figure  6  which  are  associated  with  winds 
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decreasing  speed  with  increasing  height.   The  curves  in  figu, 
5  and  6  are  similar,  except  that  persistence  occurs  only  on 
the  positive  error  side  (because  winds  are  decreasing). 
Figure  6  shows  somewhat  more  irregularity  than  figure  5  because 
of  the  smaller  sample  in  the  former.   Presuming  the  plus  skew 
of  the  geostrophic  thermal  wind  extrapolation  curve  is  valid, 
it  indicates  the  non-applicability  of  the  simple  geostrophic 

rmal  wind  for  large  curvature  cases  -  a  not  unexpected 
result  from  theory  in  view  of  the  high  percentages  of  cyclonic 

es  considered  here. 
Figure  7 

Figure  7  depicts  the  frequency  of  speed  errors  for 
cases  of  winds  increasing  with  increasing  height.   Here  is 
shown  a  pronounced  skew,  with  a  mode  at  -10  percent,  for  all 
techniques  except  persistence.   The  number  of  large  curvature 
cases  is  again  too  small  to  be  of  great  significance,  although 
the  curves  follow  those  of  figure  5  quite  closely. 

The  overall  skew  indicates  that  none  of  the  extrapola- 
tion techniques  increases  the  wind  speed  adequate ly  through 
the  layer,  the  computations  being  on  the  average   ten 
percent  too  low. 
Figure  8 

This  picture  indicates  the  extreme  variability  in  the 
accuracy  of  wind  speeds  extrapolated  through  a  layer  averag- 
ing 6500  feet.  Thermal  wind  computations  show  major  peaks 
at  both  zero  and  -20  percent  error  so  that  a  further  error 
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investigation  was  made  lecessary  (figures  9  and  10).   Persistence 

shows  a  broader  curve,  peaking  at  +10  percent  and  is  apparently 
of  less  significance.   Shear  computations  show  a  rather  irregular 
distribution,  with  the  mean  approaching  -20  percent  error. 

Large  curvature  cases  have  a  mode  at  +10  percent  with  a 
secondary  peak  at  -20  percent  error. 

Figure  9 

Considering  only  the  cases  of  wind  which  decreases  speed 
with  increasing  height  from  figure  8,  a  much  clearer  mode  at 
+10  percent  error  is  obtained.  These  curves  incjic^te  that  the 
computations  either  underestimate  the  decrease  of  wind  with 
height  or  increase  the  wind  on  the  average,  by  about  ten 
percent. 

The  data  for  the  secondary  maximum  at  -20  percent  error 
was  examined,  and  in  addition  to  several  low-wind-speed  and 
cross-isobar  cases,  there  are  four  that  cross  the  tropopause 
and  six  that  pass  through  the  layer  of  maximum  wind.   Some  of 
the  explanations  for  the  occurrence  of  this  type  of  error  would 
be  a  secondary  maximum  above  the  primary  maximum  in  the  actual 
wind  profile,  an  excessive  thermal  gradient  at  the  level  where 
the  extrapolation  began,  or  a  thermal  gradient  associated  with 
a  relative  minimum  just  above  the  original  level. 

The  error  curves  indicate  that  the  best  results  could  be 
obtained  by  using  persistence,  if  the  wind  was  known  to  be 
decreasing.   A  wind  at  the  top  of  the  layer  correct  within 
+  20  percent  would  be  obtained  60  percent  of  the  time  by  sub- 
tracting 15  percent  from  the  wind  at  the  known  level.   The 
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extreme  errors  (  >  50  percent)  that  may  arise,  are  in  part  attri' 
butable  to  the  large  number  of  cross-isobar  and  low-wind-speed 
cases  there. 

The  curve  for  large  curvature  cases  follows  the  parent 
curve  in  figure  8  except  for  the  indication  that  here  the  geo- 
strophic  thermal  wind  extrapolations  do  not  decrease  the  wind 
adequately. 
Figure  10. 

For  the  cases  of  winds  increasing  with  increasing  height 
in  figure  8,  shear  extrapolations  produce  a  maximum  skew  with 
a  mode  of  -30  percent  error.   Both  thermal  wind  computations 
have  modes  of  zero,  with  60  percent  of  the  geostrophic  and  66 
percent  of  the  gradient  thermal  wind  computations  lying  in  the 
range  -5  -  20  percent. 

. sistence  would  be  better,  however,  when  the  wind 
speed  is  known  to  be  increasing,  as  78  percent  of  the  winds 
lie  within  -  20  percent  of  a  meai   of  -15  percent. 

Large-curvature  cases  do  not  show  any  clear  indications  of 
relevant  trends. 
Table  1 

The  results  in  figures  5-10  indicate  that  it  is 
important  to  know  v/hether  the  winds  are  increasing  or  decreas- 
ing through  a  layer.   Table  1  was  assembled  to  predict  the 
accuracy  with  which  each  extrapolation  technique  could  give 
this  information.   Except  for  the  shear  computation  through 
the  thicker  layer,  the  correct  result  can  be  expected  60  per- 
cent of  the  time. 

26 


O    r:Cft$cS\-e-   2  0  ,'<A/o  rj        Q 


C  ft  s  i.  ' '  t   c  '?  a  s .»—  //■  a  /?  /•>  /*_ 


'<?*>-> 


jT/7£  i.  <P 


V  •-     A  V.,  ,-    OT«^«-fr/^  |X   -     ^  Vrf   /.xva'/i^M/,,J 


FIGURE  9 «     Frequency  Distributioi    o.C  v..';,'    Speed   Errors  (in  percent) 

t-fflQ  -'I'cjt  average   thickness  -   .-.'.J    _"c:-'    ___^     ■• '  "'    '     '  fat 

27 


o 


i  B 


•"  /    ■■-  '£        £T/e  a:  cj.-.  -  (  P/ .'/..  ce  as  y~\ 


^  *     ^    ^  *      £  1<TA.A.  'our   --> 


C7     •        <A    l<;.  ,,      *    v  p- -,-    ,-,  /   ^^r-Jx     ,      A  V/i 


;./ 


£-  V  --*/2  w/'V. /  ,}  v /..-.,. 


FIGURE  10.     Frequency  Distribution  of  Wind  Speed   Errors  (in  percent) 
6£00  feet  aver;  ~e  thickness  -  vr.'nd  incre;  :  ins  -..:1th  height 


28 


ALGEBRAIC  SIGN  OF 

ALGEBRAIC  SIGN  OF 

TYPTC  nw 

SPEED  C 

HANGE  CORRECT 

SPEED  CHANGE  II  G0;tRECT 

AZ 
210 

EXTRAPOLATL  1 

AV  f  0 

AV  =  0 

TOTAL 

AV  f  0 

AV  =  0 

TOTAL 

0  Ft 

21 

Geostrophic 
Thermal   Wind 

Ui 

18 

62 

17 

38 

Gradient 
Thermal  Wind 

U2 

20 

62 

20 

1V> 

38 

Shear 

ii2 

19     . 

62 

lii 

2)4 

38 

6^00  ff-H 

Geostrophic 
Thermal   »Vind 

$8 

3 

61 

20 

19 

39 

Gradient 
Thermal  iVind 

60 

3 

63 

19 

18 

37 

Shear 

U2 

7 

U9 

16 

35 

51 

TABLE  1.     Eumber  of  Cases  (percent)   of  Correct  Algebraic 
Sign  of  Speed  ExT.rapoj_aL±oi)t3  as  a  function  of   Type  of 
Extrapolation 
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6.   Conclusions 

From  this  study,  it  is  apparent  that  more  accurate  extra- 
polations in  the  vicinity  of  200  mb  can  be  made  if  the  influence 

of  various  restrictive  conditions  is  also  considered.   A 

it 
qualitative  estimate  of  the  likelihood  of  encountering  large 

errors  is  provided  by  the  proportionately  large  number  of 

cross-isobar  and  low-wind-speed  (and  possibly,  large-*curvature ) 

cases  occurring  near  the  limits  of  the  error  distribution 

graphs.   In  addition,  the  type  of  extrapolation  produces  an 

intrinsic  error.   Even  more  important   is  the  determination 

of  whether  the  wind  speed  is  increasing  or  decreasing 

with  height T 

Considering  wind  direction  alone,  for  a  2100-foot  layer, 
persistence  is  observed  to  be  as  good  as  either  of  the  thermal 
wind  extrapolations.   For  a  6500-foot  layer   the  gradient 
thermal  wind  computation  is  superior,  with  73  percent  falling 
within  -  15  degrees  and  88  percent  within  -  25  degrees  of 
actual  turning  through  the  layer.   For  large  curvature  cases 
in  a  6500-foot  layer,  both  thermal  wind  computations  produce 
a  small  error  that  averages  five  degrees  too  little  turning. 
Where  computing  time  is  limited,  persistence,  for  these  same 
cases,  will  produce  an  error  less  than  25  degrees  76  percent 
of  the  time. 

Better  results  are  obtained  for  vertical  extrapolations 

of  the  speed  using  thermal  wind  extrapolations  (thermal 

+  + 

gradient  accurate  to  -  15  percent  and  -  25  percent  75  and 
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90      .t  of  tiie  Li      pectively)  compared  to  persistence  or 
sheer  computations  for  a  2100-foot  layer.   For  decreasing  wind 
sp     Lth  height,  geostrophic  thermal  wind  computations  do  not 
decrease  the  wind  sufficiently,  whereas  for  increasing  speed 
h  height,  all  computations  average  ten  percent  low. 

Extrapolation  results  for  speed  are  least  accurate  for  a 
6500-foot  layer.   Generally,  best  results  are  obtained  for 
shear  computations  with  the  speed  error  in  the  range  -10  -  25 
percent,  70  percent  of  the  time.   For  decreasing  wind  speed 
with  increasing  height,  all  curves  peak  near  plus  ten  percent 
error  with  a  secondary  mode  near  -20  percent.   Consequently, 
persistence  minus  13  percent  is  correct  within  -  20  percent 
of  the  true  wind,  60  percent  of  the  time.   Even  so,  errors 
greater  than  50  percent  occur  20  percent  of  the  time,   Consider- 
ing only  increasing  winds  with  height,  thermal  wind  computa- 
tions peak  at  zero.   Again,  persistence  plus  15  percent  yields 
78  percent  of  the  cases  within  -  20  percent. 

thermal  wind  extrapolations  predict  the  correct 
direction  of  shear (plus  or  minus)  in  the  vicinity  of  200  mb 
60  percent  of  time  as  compared  to  50  percent  for  shear-wind 
computations . 

Other  approaches  that  are  recommended  for  further  study, 
to  obtain  more  accurate  wind  extrapolations,  are  listed  below. 
(1)   Instead  of  considering  speed  and  direction  separately, 
consider  the  magnitude  (and  perhaps  direction)  of  the  vector 
error  (especially  for  low  wind  speeds). 
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. 
of  large  curvature  cas.  . 
(3;        cross-isobar  cases  separately, 
(4)  re  accurate  method  of  predicting  whet 

the  wind  is  increasing  or  decreasing  with  increasing  height. 
(  i)   Separately  determine  the  shear  above  and  below  the  LMW 
and  use  for  extrapolation  studies. 
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